Kelvin-Tkachenko waves of few-vortex arrays in trapped Bose-Einstein condensates 
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We have calculated the low-lying elementary excitations of three-dimensional few-vortex arrays 
in trapped Bose-Einstein condensates. The number of different Kelvin-Tkachenko vortex wave 
branches found matches the number of vortices in the condensate. The lowest odd-parity modes 
exhibit supcrfluid gyroscopic vortex motion. Experimentally, these modes could be excited and 
observed individually or in connection with the formation and decay of quantum turbulence. 
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I. INTRODUCTION 

Rotation exerted on a superfluid forces topological 
transmutations in its structure due to the nucleation 
of quantized vortices pQ. Such remarkable behavior 
sets quantum liquids including superconductors, super- 
fluid helium and gaseous Bose-Einstein condensates apart 
from Newtonian fluids. Vortex degrees of freedom intro- 
duce additional normal modes in the spectra of superfluid 
systems altering their dynamical properties. Excitations 
of multi- vortex systems are pivotal to the understanding 
of the microscopic mechanisms which govern quantum 
turbulence in superfluid systems [21 [3] . Vortex waves are 
also relevant to the dynamics of quenched quantum phase 
transitions and such excitations are likely to accompany 
spontaneously formed vortices in the birth of superflu- 
ids |U [5]. Furthermore, the elementary excitations of 
three-dimensional rotating vortex lattices may play role 
in post-glitch relaxation phenomena observed in rotating 
neutron stars whose inner parts are thought to be filled 
by neutron superfluids [£]. 

Quantized vortex lines in simple rotating superfluids 
arrange themselves into the same triangular lattice struc- 
ture as predicted for superconductors by Abrikosov [7]. 
Kelvin waves are helical modes which propagate along 
the vortex line. They have been extensively studied in 
classical fluids since the seminal work by Thomson (Lord 
Kelvin) 8 . Their dispersion relation has the characteris- 
tic long- wavelength form wk ~ Tfc 2 ln(l/|fc|r c ), where T 
denotes the circulation around the axis of the vortex, r c 
is a core parameter, and k is the wave vector. Similar dis- 
persion relation holds also for Kelvin waves of quantized 
vortex lines in superfluids [HI HO] . 

Vortex lattices support collective Tkachcnko waves in 
which the individual vortices trace elliptical trajecto- 
ries around their equilibrium positions in the plane per- 
pendicular to the axis oriented along the vortex lines 
[TT1 IT2] . In contrast to axial Kelvin waves, transverse 
Tkachenko vortex waves have thus far only been ob- 
served in supcrfluid systems. The generic dispersion re- 
lation for Tkachenko waves in uniform systems, 

C- C^prC 



--■- '^/(clk 2 + 4fi 2 ), where c s and ct, respectively, de- 
note the propagation speed of sound and Tkachenko 
wave, changes its character from linear (stiff limit) to 



quadratic (soft limit) as the rotation frequency £1 of the 
system is increased [13] . Tkachenko waves in helium su- 
perfluids were discussed using continuum theories by Fet- 
ter and Williams [HI [TS] , Sonin [T31 HH] and Baym and 
Chandler [T7] . The stiff modes were also detected in ex- 
periments [18 . 

Observations of vortices and their dynamical proper- 
ties become particularly transparent in weakly interact- 
ing quantum gases. Orbital motion of a single quan- 
tum vortex was imaged in a Bose-Einstein condensate of 
sodium gas by Anderson et al. [19 . Such vortex mo- 
tion is a manifestation of the fundamental Kelvin wave. 
Kelvons of higher axial wave numbers were created and 
observed by Bretin et al. in elongated Bose-Einstein con- 
densates [2n|. Further discussion of those experiments 
can be found e.g. in Refs. [TUl l2"Trf2"3"] . Direct imaging of 
planar Tkachenko waves in rotating Bose-Einstein con- 
densates were reported by Coddington et al. [24] ■ Baym 
explained the observed excitation frequencies by calculat- 
ing the Tkachenko modes in two-dimensions using contin- 
uum elasticity theory accounting for the compressibility 
of the condensate [25], see also Refs. [26H30] , 

Here we obtain, for the first time, a complete mi- 
croscopic description of the lowest excitation modes 
of few-vortex arrays in harmonically trapped rotating 
Bose-Einstein condensates by employing the Bogoliubov- 
de Gennes wavefunction formalism in fully three- 
dimensional configurations. We explicitly account for 
both Kelvin and Tkachenko contributions to the collec- 
tive vortex motion and mode frequencies and include 
finite-size effects such as harmonic confinement and bend- 
ing of vortex lines. Together with sound waves, such 
elementary excitations exhaust the low-lying collective 
modes in a vortex lattice system. They are also the key 
to understanding more complex phenomena such as vor- 
tex interconnections inherent in quantum turbulence. 



II. MODEL 

The zero temperature Bose-Einstein condensate 
ground state wavefunction <p(r) and its dynamics is mod- 
eled remarkably well by the Gross-Pitaevskii equation 



ihd t </)(r,t) = C(r,t)(j)(r,t) 



(1) 



where 



C(r,t) = - 



ft 2 v s 

2m 



+ V ext (r)+ 5 |<Xr)| 2 -SlL 2 



(2) 



In Eq.|2| the constant g determines the strength of 
s— wave interactions between particles of mass to, L z 
is the component of the angular momentum operator 
parallel to the rotation axis and the external potential 
Ksxt(r) = m[w\(x 2 + y 2 ) + uj 2 z 2 ]/2 is fixed by the trans- 
verse wj_ and axial uj z harmonic frequencies. The norm 
of the wavefunction, J v \(f>(r)\ 2 dr = N, yields the num- 
ber of particles N in the system. We elect fku± and 
ao = ^/h/mu!± to be the units of energy and length, 
respectively. Our computational system is thus spec- 
ified by the dimensionless parameters uj z — 0.2 wj_, 
gN = 1000 hLu±a,Q, and the trap rotation frequency fi. 

The collective modes of the condensate are solutions 
to the microscopic Bogoliubov-de Gennes equations 

[£(r) - n]u q (r) + M(r)v q (r) = E q u q {r) 
[£*(r)-n]v g (v)+M*(r)u q (r) = -E q v q (r) (3) 

where \x = (£(r)), A4(r) — (f> 2 (r) and u q (r) and V q (r) 
are the usual quasiparticle amplitudes and E q = huj q 
the corresponding energies, where q uniquely labels the 
eigenstates. We discretize Eqs.pl) using a finite-element 
discrete variable representation [3T1 I32j and diagonalize 
the resulting eigenproblem using a parallelized implemen- 
tation of an iterative Arnoldi process, for details see [55] , 
The collective oscillatory motion of the condensate per- 
turbed by the quasiparticle modes may be visualized by 
the time-dependent probability density |</) g (r,t)| 2 with 

0,(r,t) - 4>(r)+p[u q (r)e-^ t + v* q (r)e^ t ] (4) 

and the number p controls the amplitude of the quasipar- 
ticle excitation mode relative to the ground state [34] , 



III. ROTATING GROUND STATES 

Figures [iTa)-(h) display rotating ground states calcu- 
lated using Eq.pJ) for different external angular rotation 
frequencies of the harmonic trap. Both top view (upper 
row) and side view (lower row) of the condensate density 
isosurfaces are shown. Values for the reduced rotation 
frequency 51' = tt/ujj_ and the angular momentum per 
particle ll = (L Z )/NH are marked in the frames. In the 
states shown, vortex bending is significant only near the 
condensate boundaries, however, for L' < 1 the single 
vortex state is strongly bent. For our parameters, the 
five vortex ground state assumes the shape of Olympic 
rings instead of a five-fold symmetric ring. Six vortex 
array has a single central vortex and five satellite vor- 
tices and seven- vortex state, Fig.HTe), completes the full 
two "orbital" lattice. With increasing values of W the 
centrifugal effect causes the cloud to expand in the ra- 
dial direction accompanied with the corresponding ax- 
ial shortening such that for ft' — 0.98 the shape of the 



cloud becomes roughly spherical, see Fig. fijh). Eventu- 
ally at high enough rotation frequencies the cloud enters 
the two-dimensional regime where lowest Landau level 
physics becomes relevant. 



IV. COLLECTIVE MODES 

When the system is pierced by quantized vortices each 
added vortex degree of freedom introduces new collec- 
tive modes in the system. The lowest Kelvin- Tkachenko 
mode reduces to the fundamental Kelvin mode in the 
case of a single vortex in the slow rotation limit and 
to the planar Tkachenko mode in the rapidly rotating 
(two-dimensional) system. Figure El diagrammatizes the 
different vortex wave branches for few- vortex arrays com- 
prising of one through four vortices. Diagrams on rows 
(l)-(4) apply to the respective ground states shown in 
Fig. [ija)-(d) and the columns correspond to different 
Kelvin- Tkachenko mode branches, the number of which 
matches the number of vortices in the system. The bul- 
lets mark the relative orientation of the vortices with 
respect to their equilibrium position such that the num- 
bers inside them enumerate the vortices in the counter- 
clockwise order. The two insets in Fig. [2] exemplify the 
interpretation of the diagrams by displaying the loca- 
tions and directions of motion of the vortices in the x — y 
plane. Different axial wave numbers n = 0,1,2... are 
implicitly included in each diagram. All modes in the 
first column describe in-phase oscillations of the vortex 
array as a whole while the last diagram of each row cor- 
responds to the torsional out-of-phase Kelvin- Tkachenko 
modes. The diagrams in Fig. [2] represent a qualitatively 
correct idealization of the actual vortex motion in which 
the vortices rotate along non-circular trajectories in the 
direction indicated by the arrows |34j . Denoting the dis- 
placement phase difference between any two neighboring 
vortices by Aip, the N v distinct vortex wave branches sat- 
isfy Aip = {0, 1, • • • , N v — l}|f-j where N v is the number 
of vortices in the single orbital array. 



A. Single vortex 

In the univortex case the sole vortex wave branch cor- 
responds to the Kelvin waves of the quantized vortex 
|10) . The lowest even parity n — quasiparticle exci- 
tation is the fundamental Kelvin mode which manifests 
itself as an orbital motion (frequently referred to as vor- 
tex precession in the literature) of the vortex line around 
its equilibrium position [T^l f3"5H3"?] . The lowest odd par- 
ity n = 1 Kelvin mode corresponds to the gyroscopic 
precession of a tilted vortex with a 7r phase difference be- 
tween the ends of the vortex line J33]. Such gyroscopic 
vortex dynamics was experimentally observed by Haljan 
et al. [38] and Hodby et al. [39], see also discussion 
by Stringari |40| . Modes with higher wave numbers n 
trace the rest of the Kelvin wave spectrum with increas- 
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FIG. 1. (Color online) Rotating ground states for different external angular rotation frequencies Q' . Frames in the upper 
row show top x-y view of the condensate density along the rotation axis while the frames in the bottom row display the 
corresponding side x-z view perpendicular to the rotation axis. The frame sides are 14 ao and 24 ao in extent. 



B. Double vortex 
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FIG. 2. Schematic representation of the Kelvin-Tkachenko 
mode branches. Each diagram corresponds to different vortex 
wave branch. The two insets clarify the physics encoded in 
the diagrams. Vortices denoted by the bullets are labelled 
by numbers in counter-clockwise order. The arrows next to 
the numbered bullets indicate the direction of local vortex 
motion. 



With two vortices in the condensate their dynamics is 
enriched by the vortex-vortex interaction. There are two 
distinct Kelvin-Tkachenko branches corresponding to the 
two diagrams shown in the row (2) in Fig. ^1 The first 
one comprises the "acoustic" modes which exhibit an in- 
phase, Aip — 0, elliptically polarized motion of the two 
vortices around their respective equilibrium positions. 
The second branch is formed by the "optical" modes in 
which the orbital motion of the vortex pair is out-of- 
phase with a phase difference Aip — n. In the latter 
radial motion of the vortices is strongly suppressed due 
to the angular momentum conservation. Consequently 
the dynamics appears as an azimuthal torsional vibra- 
tion [34| . Interestingly, for low momenta the modes in 
the optical branch have slightly lower excitation energies 
than the acoustic ones. 



C. Triple vortex 



ing frequencies |10] , An external rotation applied to the 
system stabilizes these Kelvin modes [4TJ [42] . As a re- 
sult, all excitation energies are positive with respect to 
the condensate energy and therefore such Kelvin modes 
rotate in the opposite sense with respect to the super- 
flow of the condensate. If the stabilizing rotating drive is 
stopped, some of the lowest modes are shifted to negative 
energies and therefore those modes change their sense of 
rotation. The Kelvin modes of an isolated vortex in a 
trapped Bose-Einstein condensate can be excited via res- 
onant coupling to quadrupole mode as detailed in Refs. 
HJJII2D1I2T1I231. 



In the three vortex ground state the vortices arrange 
themselves into the elementary triangular unit cell of the 
Abrikosov vortex lattice as shown in Fig. [Tic). Corre- 
sponding to the three vortex degrees of freedom, there 
are three branches of vortex waves [34 . In the n = 
mode of the first branch, Aip = 0, all three vortices travel 
in unison around their respective equilibrium positions 
roughly preserving their triangular relative orientation. 
The first odd mode n = 1 of this branch describes gy- 
roscopic vortex motion of the three vortex array as a 
whole. Such tilting modes were observed experimentally 
by Smith et. al |43j . The schematic in the upper in- 
set of Fig. [2] illustrates the motion of the three vortices 
in the second branch with Aip = 4ir/3. This kind of 
three-vortex oscillatory motion was observed by Yarm- 



chuck and Packard who produced photographic movies 
of few- vortex arrays in superfiuid helium [44 . The third 
branch is comprised of out-of-phase, Aip — 2tt/3, tor- 
sional Kelvin-Tkachcnko modes. As for the case of two 
vortices this last Kelvin-Tkachenko branch, referred to 
as twisting modes by Chevy 45i , contains the mode with 
lowest oscillation frequency. 



D. Quadruple vortex 



odd parity n = 1 states in the first column of Fig. [2] are 
plotted. The faster the cloud rotates the more energy 
(torque) is required to tilt the axis of rotation of the 
system reflecting the fundamental feature (conservation 
of orbital angular momentum) of gyroscopes to maintain 
their orientation with respect to the distant stars. The 
lowest axial dipole modes whose frequencies at 0.2 hu± 
are unaffected by the rotation are joined by the horizontal 
line. 



Four vortex array supports four branches of Kelvin- 
Tkachenko vortex waves [33]. There are an in-phase 
Aip = and out-of-phase A<p — it/2 branches respec- 
tively corresponding to the first and last columns of the 
row (4) in Fig. pi The motion in the second branch, 
Aip = it, may be viewed as that of a pair of vortex 
molecules such that the internal motion of vortices in 
each molecule is in-phase while the relative motion be- 
tween the two vortex molecules has a 7r phase difference. 
Vortex motion in the remaining third branch is illustrated 
in the lower inset of Fig. [2J In it Aip — 3tt/2 and the 
internal motion of each vortex molecule is out-of-phase 
and also the relative motion between these two vortex 
molecules has a 7r phase difference. 



E. Spectra 

Figure [3] (a) shows the energies of the lowest quasi- 
particle excitation modes for the three-vortex ground 
state of Fig. [T] as functions of their axial momenta de- 
fined by k z = VK0)l - d 2 z \u q {v)) + (v q (r)\ - d 2 z \v q {r)). 
The three different Kelvin-Tkachenko branches are la- 
beled using the markers indicated in the figure. The 
lowest odd-parity modes of each branch correspond to 
the tilted gyroscopic motion of the vortices. Above the 
Kelvin-Tkachcnko branches in energy lie the sound waves 
with their characteristic linear dispersion relation which 
becomes curved due to the density inhomogeneity. The 
lowest co-rotating centre-of-mass mode has the frequency 
uj± — £1 and at Q = lj±_ it grows without a bound causing 
the whole condensate to spiral out of the trap due to the 
centrifugal weakening of the effective trapping potential. 

In Fig. |3Fb) we have plotted the n — modes of each 
Kelvin-Tkachenko branch with markers ▲, I, T, and ♦ 
corresponding to the diagrams in the last, second but 
last, third but last and fourth but last column of each row 
of Fig(2] respectively. These modes tend toward the zero 
point energy as O is increased becoming the constituents 
of the lowest Landau level in the rapid rotation limit. In 
addition, the gyroscope modes corresponding to the first 



V. CONCLUSIONS 

To conclude, we have computed the low-lying elemen- 
tary excitation modes of three-dimensional few- vortex ar- 
rays in rotating trapped Bose-Einstein condensates. The 
lowest energy Kelvin-Tkachenko excitations correspond 
to an oscillating torsional motion of the vortex array in 
the plane perpendicular to the rotation axis. The lowest 
odd axial modes exhibit superfiuid gyroscopic vortex dy- 
namics. The Kelvin-Tkachenko waves could be created 
in ultra-cold atomic gases experimentally using a suit- 
able combination of laser potentials [20, 24 or they might 
become spontaneously excited in quenched phase transi- 
tions 4, 5 . Continuous imaging methods should allow 
experimental investigation of the dynamics of these vor- 
tex waves [5] . In superfiuid helium vibrating wire exper- 
iments [3S1 I3Z] might enable experimental investigation 
of such elementary excitations. 

In multi- vortex systems vortex-vortex interactions en- 
rich the internal dynamics of the superfiuid. As 
in the case of a single vortex |10j . time-dependent 
Gross-Pitaevskii simulations show that certain Kelvin- 
Tkachenko waves of multi- vortex arrays can be excited by 
resonantly populating the counter-rotating quadrupole 
mode which couples to these vortex waves. Resonant am- 
plification of multi-vortex Kelvin-Tkachenko waves can 
lead to vortex reconnections when the amplitude of the 
modes exceeds the intervortex spacing, onsetting quan- 
tum turbulence in these superfiuid systems. In future, 
systematic investigation of interaction dynamics of vor- 
tex waves in Bose-Einstein condensates may shed light on 
the microscopic mechanisms which govern quantum tur- 
bulence such as vortex reconnections, cascades and the 
role of Kelvin-Tkachenko waves in transport and dissipa- 
tion of turbulent energy. 
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